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Abstract. Giant planets found orbiting close to their central stars, the so-called 'hot Jupiters', are thought to have originally 
formed in the cooler outer regions of a protoplanetary disk and then to have migrated inward via tidal interactions with the 
nebula gas. We present the results of N-body simulations which examine the effect such gas giant planet migration has on 
the formation of terrestrial planets. The models incorporate a 0.5 Jupiter mass planet undergoing type II migration through an 
inner protoplanet-planetesimal disk, with gas drag included. Each model is initiated with the inner disk being at successively 
increased levels of maturity, so that it is undergoing either oligarchic or giant impact style growth as the gas giant migrates. 
In all cases, a large fraction of the disk mass survives the passage of the giant, either by accreting into massive terrestrial planets 
shepherded inward of the giant, or by being scattered into external orbits. Shepherding is favored in younger disks where there 
is strong dynamical friction from planetesimals and gas drag is more influential, whereas scattering dominates in more mature 
disks where dissipation is weaker. In each scenario, sufficient mass is scattered outward to provide for the eventual accretion of 
a set of terrestrial planets in external orbits, including within the system's habitable zone. This scattering, however, significantly 
reduces the density of solid material, indicating that further accretion will occur over very long time scales. A particularly 
interesting result is the generation of massive, short period, terrestrial planets from compacted material pushed ahead of the 
giant. These planets are reminiscent of the short period Neptune-mass planets discovered recently, suggesting that such 'hot 
Neptunes' could form locally as a by-product of giant planet migration. 
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1. Introduction. 

Over the last ten years the radial velocity technique has been 
successfully employed in the detection of giant planets orbit- 
ing nearby main sequence stars (e.g. Mayor & Queloz 1995; 
Butler et al. 2004; Marcy, Cochran & Mayor 2000). To date, 
136 extra-solar planetary systems have been discovered, 14 of 
them multiple, leading to a total of 155 giant planets'. Even 
though radial velocity observations are more sensitive to short 
period orbits, a surprising discovery has been the substantial 
population of giant planets orbiting close to their central star at 
distances < 0.1 AU. Twenty nine such objects are known (3 of 
them sited in multiple systems) comprising ~ 20% of the total 
sample. These so-called 'hot Jupiters' are mostly sub-jovian in 
mass, with low eccentricity orbits, and are often found associ- 
ated with stars more metal rich than the Sun (Udry et al. 2003; 
Santos et al. 2003; Fischer & Valenti 2005). 

The most likely scenario for the origin of hot Jupiters in- 
volves initial formation further out in the nebula beyond the 
'snowline' as per conventional formation theories (e.g. Pollack 
et al. 1996), followed by an episode of inward orbital migra- 
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tion, propelled by tidal interactions between the planet and gas 
disk. A variety of orbital migration phenomena have been pro- 
posed for protoplanets as they grow in mass, but the one most 
likely to operate in this case is type II migration where the 
planet has become sufficiently massive (<: 100 M®) to open up a 
gap in the gas, thereby migrating inward at a rate controlled by 
the disk viscous time scale (Lin et al. 1996; Lin & Papaloizou 
1986; Nelson et al. 2000). This time scale is typically on the 
order of a few xlO^ years. What finally stops the inward drift 
of these planets at such small radial distances, other than for- 
tuitous disk dispersal, is presently unknown and it may be that 
some migrating planets are accreted by the central star. 

It appears therefore that hot Jupiter systems are not uncom- 
mon and possibly represent an extreme rearrangement of plane- 
tary mass during formation as compared with the Solar System. 
This comparison raises an interesting set of questions. What ef- 
fect would hot Jupiter migration have on teiTestrial planet for- 
mation? Might accretion of teiTestrial planets be inteiTupted 
or prevented altogether? Can we have any realistic expecta- 
tion of ever discovering Earth-like planets in these systems, 
or is it probable their habitable zones are empty of significant 
material, having been swept clean by the passage of the gi- 
ant? The assumption that hot Jupiter systems are barren are 
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among those advanced to support such speculative concepts as 
the "Rare Earth" hypothesis (Ward & Brownlee 2000) and the 
Galactic Habitable Zone (Lineweaver 2001; Lineweaver et al. 
2004), both of which argue for a large number of special cir- 
cumstances required to form an Earth-like planet. From their 
standpoint, any significant diff'erence in planetary system ar- 
chitecture from that of the solar system is Ukely to preclude 
the formation or survival of planets that are habitable for multi- 
cellular hfe. 

These issues were examined by Armitage (2003) with a 
simple time-dependent disk model of the evolution of gas, dust 
and planetesimals. He assumed that the eff'ect of giant planet 
migration is to first sweep the inner disk clear of planetesimals, 
and then looked at whether a new generation of planetesimals 
could be formed in the terrestrial planet region by the subse- 
quent resupply of dust from outer regions by advection and dif- 
fusion. His conclusion was that resupply of solid material into 
the inner disk would be ineflicient and terrestrial planet for- 
mation would be unlikely. In contrast, Mandell & Sigurdsson 
(2003) consider a late migration scenario and use N-body simu- 
lations to model the migration of a Jupiter mass planet through 
a fully formed terrestrial planet system. The typical pattern of 
evolution observed included: 1) excitation of planetary orbits 
by sweeping resonances with the inward migrating giant; 2) 
close encounters between the planets resulting in colUsions or 
mutual scattering; and 3) slingshot encounters with the giant 
as it passed through the inner system, leading to ejection, col- 
Usion with the central star, or scattering into bound exterior 
orbits with increased semi-major axis, eccentricity and inclina- 
tion. Short migration times allowed a larger fraction of planets 
to survive rather than being ejected whereas the trend was re- 
versed for long migration times. Overall, ~ 25% of the planets 
survived in a wide variety of orbits exterior to the giant, orbits 
which, they speculated, might subsequently become circular- 
ized as a result of dynamical friction with outer system plan- 
etesimals or interaction with the remnant gas disk. They con- 
cluded therefore that inward migration of a giant planet does 
not invariably eliminate pre-formed terrestrial planets and that, 
given an initial layout of bodies similar to that of the Solar 
System, between ~ 1-4% of systems in which migration oc- 
curred could still possess a planet in the habitable zone. Actual 
formation of terrestrial planets in the presence of a hot Jupiter 
has been modelled by Raymond et al. (2004). They assume a 
previous rapid migration of the giant into its final close orbit 
and then model the later stages of terrestrial planet accretion 
from an exterior protoplanet disk using N-body methods. Their 
conclusion is that the presence of a hot Jupiter does little to in- 
terfere with terrestrial planet formation outside of an annulus 
that is within a factor of three in period to the giant (about a 
factor of two in semi-major axis). Planet formation in the hab- 
itable zone is not adversely affected. 

The conclusions of these three papers span the widest pos- 
sible range of outcome, from the occurrance of terrestrial plan- 
ets in hot Jupiter systems being highly unlikely, through pos- 
sible but rare, to commonplace. This variation originates from 
the very different assumptions and initial conditions in each 
case. Armitage (2003) does not model the dynamic effects of 
the migrating giant on planetesimals and instead assumes to- 



tal loss of planetary building blocks from within the swept 
zone. This may be unrealistic especially as, by the time a giant 
planet has grown large enough to start type II migration, con- 
siderable accretion into larger planetary embryos could have 
already occurred in the inner system, bodies which might not 
be so readily accreted or bulldozed into the central star. The 
picture of Mandell & Sigurdsson (2003) of a giant migrating 
through a mature terrestrial planet system may suffer from un- 
realistic timing as giant planet formation and migration is con- 
strained to occur within the ~ 10^-10^ year lifetime of the gas 
disk whereas the terminal 'giant impacts' phase of terrestrial 
planet formation is thought to last ~ 10^ years (e.g. Chambers 
2001). The optimistic conclusions of Raymond et al. (2004) are 
a function of their initial condition of placing the hot Jupiter 
in its final, post migration, orbit and then modelling terrestrial 
planet formation in an essentially undisturbed, unmixed, exte- 
rior disk. This requires rapid giant planet formation and migra- 
tion followed by the formation of a new terrestrial disk from 
the remaining debris which must somehow be reduced back to 
a more unevolved, damped, and chemically differentiated state. 

Timing issues are therefore important in the study of this 
problem, involving initially isolated sequences of events inside 
and beyond the nebula snowline. Whilst there remains much 
uncertainty over detail, one fairly certain constraint which gives 
an upper Umit to the time available is that giant planets must 
both form and complete migration in considerably less than 10^ 
years, before the dispersal of the nebular gas. Observations sug- 
gest that 50% of stars in clusters have lost their disks by an age 
of ~ 3 X 10^ years (Haisch et al. 2001). Estimation of a lower 
age limit is more problematic as it must rely on our incom- 
plete theories of giant planet formation. If the core-accretion 
model is to be preferred, favorable conditions would allow gi- 
ants to form in > 10^ years (Pollack et al. 1996; Papaloizou & 
Nelson 2004; Alibert et al. 2004). In the meantime, accretion 
will be ongoing within the planetesimal swarm in the terres- 
trial planet region. According to the current picture, an early 
phase of runaway growth will give way to a more lengthy phase 
of oligarchic growth where similar sized protoplanets emerge 
from the swarm in well-spaced orbits which remain near circu- 
lar due to dynamical friction from the surrounding sea of plan- 
etesimals (Kokubo & Ida 1998). Oligarchic growth ends when 
planetesimal numbers decline to the extent that their damping 
effect on protoplanet orbits becomes insufficient to prevent or- 
bit crossing. This inaugurates the last phase of terrestrial planet 
formation, that of so-called 'giant impacts', involving the mu- 
tual accretion of protoplanets and thinning down of their num- 
ber to the point where the final planets emerge, positioned in 
stable non-crossing orbits. Simulations of this final stage of ter- 
restrial planet growth suggest that it would take ~ 10^ years 
to complete (Chambers 2001), long after the disappearance of 
the nebular gas. Oligarchic growth however starts much ear- 
Uer, whilst gas is still present: simulations by Kokubo & Ida 
(2000) have shown that it takes only ~ 5x10^ years to gen- 
erate ~ 0.01 - 0.03 Me planetary embryos from a planetesi- 
mal disk at 1 AU. Thus, in the case of a giant planet migrating 
through the terrestrial planet zone, it seems most probable that 
this would occur at some time within, or towards the end of the 
phase of oUgarchic growth in that region. 
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In order to better satisfy these timing constraints, 
we present below a set of N-body simulations of giant 
planet migration through progressively evolved inner sys- 
tem protoplanet-planetesimal disks undergoing either late oli- 
garchic or early giant impact style growth. In Section 2 we out- 
Une our model and its initial conditions; in Section 3 the results 
are presented and discussed; in Section 4 we consider some 
caveats and future model improvements, and in Section 5 we 
offer our conclusions. 



2. Description of the model. 

2.1. Forces. 

We choose to model planet growth and migration using 
the hybrid- symplectic N-body simulation package Mercury 6 
(Chambers 1999), modified to include the effect of gas drag on 
planetesimals and type II migration on a single giant planet. 
The ingredients required for the simulations are thus one mi- 
grating giant planet, and an interior disk of small protoplanets 
embedded within a swarm of planetesimals. However, due to 
the huge number of particles involved, the realization of a re- 
alistic planetesimal disk, with every body treated as fully in- 
teracting, is well beyond the current state of the art. We pro- 
ceed therefore by following the so-called N + N' approach of 
Thommes et al. (2003) where we have protoplanets embed- 
ded in a disk of A^' 'super-planetesimals', particles that rep- 
resent an idealized ensemble of a much larger number of real 
planetesimals. The giant and the protoplanets feel all the mod- 
eled gravitational forces, whereas the super-planetesimals feel 
gravitational forces from the central star, protoplanets, and gi- 
ant planet, but are otherwise non self-interacting. This prevents 
their relatively high masses from unrealistically auto-exciting 
the disk. Super-planetesimals alone also experience gas drag 
with the drag being calculated using a defined physically real- 
istic planetesimal radius. The issue of the mass of the super- 
planetesimals was addressed by Thommes et al. (2003) who 
found in test runs that protoplanets undergo effective dynam- 
ical friction if super-planetesimal masses are ^ 0.1 times the 
initial protoplanet masses. 

The coordinate origin is based on the central star. The accel- 
eration experienced by each of the super-planetesimals is given 
by: 
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where M, is the stellar mass, mj and rj are particle masses and 



position vectors, and r,j 



The first term on the right 



hand side represents the acceleration from the central star, the 
second term the accelerations from the protoplanet and giant 
planet, the third and fourth terms are the indirect terms arising 
from the acceleration of the coordinate system, and ttdrag is the 
acceleration due to gas drag. 

The acceleration experienced by the protoplanets and gas giant 



planet is given by: 
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The first term again represents the acceleration due to the cen- 
tral star, the second and third terms represent the accelerations 
due to the protoplanets/giant planet and super-planetesimals, 
respectively, and c5,j is the Kronecker delta function. The fourth 
and fifth terms are indirect terms, and axypeii is the acceleration 
driving type II migration of the gas giant planet alone. 

Planetesimals are small enough to experience a drag force 
from moving through the nebula gas. This acceleration, which 
acts both to cause an inward radial drift and a damping of ec- 
centricities and inclinations, takes the form: 
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where mpi and rp\ are the physical mass and radius of a single 
planetesimal respectively, Co is the drag coefficient (taken here 
to be Cd - 1) and u = Vp - Vg is the velocity of the planetesimal 
with respect to the gas. The gas is assumed to move in a circular 
orbit which, due to pressure support, revolves at slightly less 
than Keplerian speed. The relation is: 



Vg = (1 - 2r])2 vk 



(4) 



where Vk is the local Keplerian velocity and 77 - 
0.0019(a/lAU)^/^, given the nebula scale height introduced 
later in Eq. 9 (Adachi et al. 1976). 

In order to model type II migration of a giant planet, we 
adopt a simple prescription and assume that this occurs at a 
rate controlled by a local viscous disk evolution timescale that 
is proportional to orbital period. Assuming an alpha viscosity 
model, this is roughly: 



(5) 



where r is the radial distance, h is the disk scale height, a is 
the alpha parameter and Q is the orbital frequency. Given a - 
2 X 10~^ h/r = 0.05 and r = 5 AU then Ty ^ 0.25 Myr. Since 
Ty oc a^-^, where a is the semi-major axis, a oc a""-^ so inward 
migration speeds up as it proceeds. (Note that this value of h/r 
differs slightly from that obtained from Eq. 9.) 

The type 11 migration process is also assumed to exert 
strong eccentricity and inclination damping which is taken here 
to operate over a timescale = Ty/50. The acceleration axypen is 
therefore implemented as: 



«Typen 



2Ty 



-25 



2(v • r)r 2(v • k)k 



(6) 



where v is the giant planet's velocity vector and is a unit vec- 
tor in the vertical direction. Note that the factor of 2 appearing 
in the the first term on the right hand side arises because the 
migration time is half the angular momentum removal time. 
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Fig. 1. Scenario I at 20 000 years after the start of giant planet migration, showing the mass, inclination and eccentricity of objects. Small black 
dots represent super-planetesimals; white filled circles are rocky protoplanets; grey filled circles are icy protoplanets and the large highlighted 
grey filled circle is the giant. The dotted line in the upper panel shows the eccentricity at which the pericentre of an exterior object intersects 
the orbit of the giant. The location of the 2: 1, 3:2 and 4:3 resonances with the giant are indicated. 



2.2. The Nebula. 

The nebula model used for this study is based on the minimum 
mass solar nebula (MMSN) of Hayashi (1981) and is defined 
as follows. 

The surface density of solids is: 

- /neb/ice2^1 [YmJ j ' 

where /neb is a nebular mass scaling factor, 2i = 7 g cm"^ and 
the ice condensation coefficient /ice = 1 for a < 2.7 AU (the 
distance chosen for the nebula 'snowline') and /ice - 4.2 for 
a > 2.7 AU. 

The volume density of gas is: 

u 

Pg=/nebPl(^ ' exp[-z2//z2j , (8) 



where pi = 2.0 x 10-'^(/gas/240)(Ii/10) g cm-^ /g^s is the gas 
to dust ratio, z is the height from the midplane of the nebula 
and the disk scale height h is taken to be: 

= 0.045 (^f. (9) 

In the context of the core-accretion model of giant planet 
formation, a significant additional mass of solids seems to be 
required in order to form a core quick enough to initiate gas 
accretion before the loss of the gaseous component of the 
nebula (Lissauer 1987; Pollack et al. 1996; Thommes et al. 
2003). Moreover, since hot Jupiters are usually found around 
stars more metal-rich than the Sun, a greater solids content 
might reasonably be expected within their protoplanetary neb- 
ulae. Thus, for this study we assume an equivalent mass of 
3 X MMSN of solids. However, as the nebula is taken to be 
somewhat evolved at our starting point, considerable gas could 
already have been lost so we assume a smaller equivalent mass 
of 2 X MMSN of gas. The relevant parameters are therefore set 
to /leb = 3 and /gas = 160. 
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Fig. 2. Scenario I at 100000 years after the start of giant planet migration. The giant has now moved inward to 2.68 AU and has scattered a 
significant amount of mass into exterior orbits. 



Table 1. Initial disk set-up 





Rocky Zone 


Icy Zone 


Total 




0.4-2.7 AU 


2.7-4.0 AU 


0.4-4.0 AU 




9.99 Me 


14.8 Me 


24.8 Me 


'^proto 


0.025 Me 


0.1 Me 




N 


66 


9 


75 




0.0025 Me 


0.01 Me 




N' 


3336 


1392 


4278 


/proto 


0.17 


0.06 


0.1 



2.3. Initial conditions and running of the simulations. 

The interior disk of protoplanets and planetesimals is modeled 
initially from 0.4—4 AU which, given the nebula parameters al- 
ready described in Section 2.2, amounts to Mjoiid - 24.8 M® of 
solid (rocky or icy) material. Adopting the oligarchic growth 
picture advanced by Kokubo & Ida (2000), we set a nominal 
age for the disk of 0.5 Myr, by which point protoplanets in the 
inner disk may have grown to a few percent of an Earth mass 
and perhaps to greater than this beyond the snowline. Initial 
protoplanet masses of wiproto — 0.025 and /Wpi-oto — 0.1 Me 



are therefore chosen to represent bodies interior and exterior 
to the snowline respectively. The number of protoplanets 
was calculated based on an assumption that the average radial 
spacing between them is 8 mutual Hill radii. Semi-major axes 
for this number of protoplanets are then generated randomly 
with probabilities weighted in order to reproduce the disk sur- 
face density profile. Eccentricities and inclinations are random- 
ized from a Rayleigh distribution with RMS values of 0.01 
and 0.005 respectively. Additional orbital elements required 
are randomized uniformly from within their range. The total 
mass of protoplanets is then subtracted from the total mass of 
the disk annulus and the mass that remains is divided into A'' 
super-planetesimals with a mass (ms_pi) one tenth that of the 
protoplanets in their zone. The orbital elements of the super- 
planetesimals are then generated in the same manner Note that 
the super-planetesimal mass ms-pi is distinct from the physical 
planetesimal mass nipi which is used solely for the purpose of 
calculating gas drag (see Eq. 3). 

Data for this initial disk model are shown in Table 1 . The 
overall values are N -15 and A^' - 4278 and the mass fraction 
of the disk contained in protoplanets /pioto ~ 0.1. The parame- 
ter /pioto is used here as a rough measure of the evolution of the 
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Fig. 3. Scenario I at 160000 years after the start of giant planet migration. The giant has now moved inward to 0.52 AU. The scattered exterior 
disk has grown, but a substantial amount of mass in planetesimals and six rapidly accreting protoplanets remain interior to the giant. The outer 
three of these six are in first order resonances with the giant. 



disk and we take /proto = 0.5, the point where the total mass in 
protoplanets exceeds that in planetesimals, to denote the tran- 
sition between oligarchic and giant impact growth regimes. 

Mercury 6 models accretion as inelastic collisions between 
objects whose radii are calculated from an input density, val- 
ues of 1, 2 and 3 g cm"^ for giant, icy and rocky masses re- 
spectively being chosen here. All collision permutations are 
enabled, except that between super-planetesimals, and a phys- 
ically realistic planetesimal radius of rpi - 10 km is imposed 
both for calculating drag and collisions. The mass of the central 
star is taken to be = 1 Mq. 

Before the giant is added to the picture, five instances of 
this disk were allowed to evolve by being run for 0.1, 0.25, 
0.5, 1.0 and 3.0 Myr, the purpose being to provide the basis for 
five type II migration scenarios through progressively evolved 
inner system material. In these runs, the timestep for the sym- 
plectic part of Mercury's hybrid integrator was set to 8 days 
and all mass straying inward of 0.1 AU was eliminated and 
added to the mass of the central star. Data for these evolved 
disks are given in Table 2. It can be seen that the amount of 
mass lost to the central star, principally via gas drag-induced 



Table 2. Overall disk data: after 0.1-3.0 Myr of evolution 



Time (Myr) 


0.1 


0.25 


0.5 


1.0 


3.0 


Scenario ID 


I 


II 


III 


IV 


V 


M,oM (M<b) 


24.4 


23.7 


22.9 


20.6 


14.6 


mmax (Me) 


0.22 


0.33 


0.52 


0.91 


1.28 


N 


56 


51 


40 


29 


19 


N' 


3863 


3312 


2660 


1341 


499 


,/proto 


0.20 


0.24 


0.32 


0.51 


0.70 



orbital decay of planetesimals, remains modest (except in the 
case of Scenario V) but that as the disk ages the maximum pro- 
toplanet mass nimax increases and particle numbers decrease. 
As would be expected, /pioto increases with time until by 1 Myr 
/proto > 0.5 and giant impact style growth has begun. 

The five type II migration scenarios studied here are con- 
structed from these five evolved disks. A giant planet of mass 
0.5 Mj is placed into each simulation at 5 AU and allowed to 
migrate inward with Ty(5 AU) = 237 000 years (see Eq. 5). 
The giant migrates down to 0.1 AU in ^ 160000 years, the 
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Fig. 4. Detail of the interior regions of Scenario 1 at 160000 years 
after the start of giant planet migration, showing e vs a for objects 
< 2 AU. Planetary masses are indicated in units of Me. A total of 
~ 15Me of material has been pushed inward by the giant, ~ two thirds 
of it remaining in super-planetesimals, many of which do not show 
individually in this diagram as they are over-plotted in the vicinity of 
the 3:2 and 2:1 resonances. 



exact time depending on the amount of mass remaining inte- 
rior to its orbit. All simulations therefore were run for 170000 
years, with the type II migration algorithm being switched 
off once the giant reached 0.1 AU. In order to better model 
processes when the giant migrates down to small radial dis- 
tances, collision with the central star is computed when r < 
0.014 AU = 3 Rq, the approximate radius of a T-Tauri star. The 
initial timestep chosen for the symplectic integrator was 1 day; 
but as dynamical spreading and the effects of migration and 
drag drives some material into closer orbits, the timestep was 
reduced as the simulation progressed. The position of the inner 
edge of the swarm was therefore monitored at intervals during 
each run keeping the timestep close to one tenth the orbital 
period of the innermost object. This considerably increased 
the run times of these simulations, especially those based on 
younger disks, requiring about a month of processor time for 
completion. 

3. Results. 

3.1. Typical features of a run. 

The character of the planetary systems formed from these runs 
was found to vary systematically with the age of the inner disk. 
However, all scenarios also exhibited a number of behavioral 
features in common. We discuss these first by describing one 
of the scenarios in detail. 

Five snapshots of the evolution of Scenario I are illustrated 
in Figs. 1-5 showing the mass, inclination and eccentricity of 
objects vs. semi-major axis. The original provenance of the 
protoplanets (interior or exterior to the snowline) is denoted by 
the shading of its symbol as described in the caption to Fig. 1 . 



In the case of a merger between rocky and icy protoplanets, this 
shading is determined by that of the most massive of the pair. 

A juncture early in the simulation is illustrated in Fig. 1, 
20000 years after the start of migration when the giant has 
moved inward to 4.58 AU. Several prominent features have de- 
veloped. The outer edge of the disk has been pushed inward at 
the 4:3 resonance with the giant - a shepherding process that 
acts primarily on planetesimals that are damped by gas drag 
(Tanaka & Ida 1999). In addition, sweeping resonances have 
captured a substantial population of planetesimals and some 
protoplanets at 3:2 and 2:1, increasing their concentration at 
these locations and exciting their orbits. The overall effect is a 
compaction and excitation of the outermost annulus of the disk 
between the 4:3 and 2:1 resonances. Inwards of this zone, the 
presence of the approaching giant has had little influence. 

The state of play some time later at 100000 years is shown 
in Fig. 2. The giant has now moved inward to 2.68 AU, con- 
tinuing to push the outer edge of the disk ahead of it at the 4:3 
resonance. An increased amount of mass has been entrained 
in the region between 4:3 and 2:1 and resonant pumping and 
mutual scattering has raised the eccentricities of some proto- 
planets to high values. This has allowed some objects to cross 
the gap between the disk and giant whereupon a close slingshot 
encounter causes expulsion into an exterior orbit. A diffuse and 
excited exterior disk is now in the process of formation, com- 
posed predominantly of the more weakly damped protoplanet 
material. 

At 160000 years the giant has moved inward to 0.52 AU 
as shown in Fig. 3. The most prominent feature now appears to 
be the scattered exterior disk made up of numerous protoplan- 
ets with large a and e and a diffuse population of planetesimals 
with a Es of only a few percent of the original disk. This how- 
ever still represents the minority of solids mass. Two thirds of 
the disk mass remains interior to the giant in the form of re- 
maining planetesimals and six protoplanets. A blow-up of the 
interior regions of the system at this juncture is shown in Fig. 4. 
A total of ~ 15Me of solid material has been compacted inte- 
rior to the giant, most of which lies between 0.32-0.39 AU. 
Two massive protoplanets of 3.57 and 1.14 are captured at 
the 3:2 resonance with the giant and a smaller 0.48 M© proto- 
planet is found at the 2:1 resonance. The majority of the mass 
however remains in super-planetesimals entrained at these res- 
onances (and mostly over-plotted in the figure) and in a ring 
of matter between them. Protoplanetary eccentricities therefore 
remain low, even at resonances, due to strong dynamical fric- 
tion. In addition, accretion rates onto these objects have now 
become so high that collisional damping is also acting to con- 
trol the growth of their eccentricities. 

The final snapshot of Scenario I at 170000 years is shown 
in Fig. 5. The giant has stopped its migration at 0.1 AU and, 
whilst the character of the exterior disk is unchanged, evolu- 
tion has proceeded rapidly to /proto = 1 within the compacted 
interior material. No planetesimals remain, having either been 
swept up by protoplanets or accreted by the central star. An in- 
tense episode of runaway accretion and giant impacts has ended 
in the assembly of a single 15.65 M© planet out of resonance 
with the giant at 0.055 AU. This arrangement of close orbiting 
giant and an inner Neptune-mass planet has a striking similarity 
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Fig. 5. Scenario I at 170000 years after the start of giant planet migration. The giant has now stopped at 0.1 AU. All interior mass has accreted 
into a single ~ 16 object after an intense episode of runaway accretion and giant impacts. 



to the two innermost planets in the 55 Cancri system (Mc Arthur 
et al. 2004). 

To reinforce the above interpretation of processes at work 
in Scenario I, the surface density evolution of the disk and its 
accretion rate are shown in Fig. 6. The left hand panel shows 
the disk surface density profile (obtained by summing all pro- 
toplanets and super-planetesimals in 0.1 AU width bins) at 
20000, 60000, 120000 and 160000 years after the start of 
migration; the right hand panel plots the amount of mass ac- 
creted onto protoplanets only (including giant impacts) every 
1 000 years for the duration of the run. In the Es plot, two sur- 
face density enhancements are clearly visible as spikes at the 
3:2 and 2:1 resonances and are seen to grow whilst moving in- 
ward. At 120000 years, these have merged into one: the outer 
half of the original disk having by now been squeezed into a 
dense ring. By 160000 years, most of this mass is now con- 
fined within 0.5 AU and Ss here has risen to ~ 10^ g cm"^ 
(an increase by a factor of ^ 10 over the previous, undisturbed, 
disk surface density) which is off the vertical scale in the fig- 
ure. The effect of this disk compaction process is seen clearly 
in the accretion rate plot. Mass accretion rises significantly af- 
ter 120 000 years due to both the resultant high values of and 



the fact that much of this mass now resides in a zone where dy- 
namical times are shorter. Growth interior to the giant ends in a 
short-lived and dramatic phase of runaway planetesimal accu- 
mulation and giant impacts within the material pushed into the 
small volume inside 0. 1 AU from the central star. 

The behavioral features seen to a greater or lesser extent in 
all the runs summarize as follows. 

1. Shepherding of planetesimals. As planetesimal random 
velocities are continuously damped by gas drag, their 
tendency is to be pushed inward, ahead of the giant. 
Shepherding of protoplanets also results as a weaker sec- 
ondary effect as they are, to a varying extent, coupled to the 
planetesimal disk by dynamical friction. 

2. Resonant capture. First order resonances with the giant 
gather an increasing amount of mass as they sweep inward. 
This, in addition to the shepherding effect eventually re- 
sults in the compacting of some of the disk mass into a 
zone close to the central star. 

3. Acceleration of planetary growth interior to the giant. 
Accretion speeds up within the compacted interior disk. 
This is particularly rapid within the disk remnant squeezed 
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Fig. 6. Surface density evolution (left hand panel) and accretion rates (right hand panel) for Scenario I. Growing surface density peaks at the 
2:1 and 3:2 resonances sweep through the inner system ahead of the giant. Accretion rates increase after ~ 120000 years and the final intense 
accretion spike represents clear up of remaining material shepherded within 0.1 AU. 



inside 0. 1 AU, where the final evolutionary phases of accre- 
tion are rushed through in mere thousands rather than mil- 
lions of years. Typically, 1-3 massive close orbiting planets 
are the end result. Where there is one survivor, its mass and 
configuration can be reminiscent of the 'hot Neptune' type 
of planet identified recently. 
4. Creation of a scattered exterior disk. Pumping of eccentric- 
ities at resonances and by mutual perturbations permit some 
disk material to undergo a close encounter with the giant 
where it is scattered into an external orbit. Protoplanets, 
being less strongly damped, are more likely to have this 
happen than planetesimals. The result is a dynamically ex- 
cited and widely dispersed external disk of material where 
accretion rates are greatly reduced. Leaving aside the issue 
of interaction of this ejecta with the disk of gas and solids 
outside the giant's original formation orbit (not simulated 
here), it seems likely from the result shown in Fig. 5 that 
further giant impact style evolution ending with a set of 
planets in non-crossing orbits would take much longer than 
the ~ 10** years estimated for the solar system. 



3.2. Dependence on the maturity of the inner disk. 

The purpose of running five scenarios through a progressively 
more mature inner disk is to explore the issue of whether the 
timing of migration has any systematic effect on the results. 
This is possible because when a disk evolves and small objects 
accumulate onto larger ones, both dynamical friction and gas 
drag become less effective overall, influencing both the shep- 
herding and scattering behaviors previously described. 

In all the simulations there were five possible fates await- 
ing all the modeled disk particles: 1) survival in a body orbiting 
interior to the giant; 2) survival in a body orbiting exterior to 
the giant; 3) accretion by the central star; 4) accretion by the 
giant; and 5) ejection from the system. These data are shown in 
Table 3 which lists the fate of the disk mass at the end of the 
simulations: the total surviving solids and the five end points 
being shown as a percentage of the total initial solids. It is no- 
ticeable from Table 3 that systematic trends in the fractionation 
of solids between some of these end points do appear 

The most obvious trends are that the percentage of the orig- 
inal disk mass that ends up surviving interior to the giant falls 
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Table 3. Fate of the disk mass at 170 000 years 



Scenario 


1 


11 


111 


IV 


V 


Total Initial Solids (Me) 


24.39 


23.68 


22.89 


20.61 


14.59 


Total Surviving Solids (Me) 


22.06 (90%) 


21.04(89%) 


12.05 (53%) 


14.75 (70%) 


9.91 (68%) 


Interior Surviving Solids (Me) 


15.65 (64%) 


14.60 (62%) 


4.96 (22%) 


2.86 (14%) 


0.00 (0%) 


./proto 


1, 1 


3, 1 


1, 0.99 


1, 1 


0,0 


Exterior Surviving Solids (Me) 


6.41 (26%) 


6.44 (27%) 


7.09 (31%) 


11.59 (56%) 


9.91 (68%) 


yproto 


32, 0.66 


32, 0.63 


25, 0.65 


23, 0.93 


16, 0.89 


Accreted by Star (Me) 


1.05 (4%) 


0.74 (3%) 


9.60 (42%) 


5.81 (28%) 


3.03 (21%) 


Accreted by Giant (Me) 


1.19(5%) 


1.88 (8%) 


1.24 (5%) 


0.34 (2%) 


1.65 (11%) 


Ejected (Me) 


0.00 (0%) 


0.025 (0.1%) 


0.00 (0%) 


0.025 (0.1%) 


0.00 (0%) 



100 




Fig. 7. Interior and exterior surviving solids as a percentage of initial 
disk mass at 170000 years. 



orbital decay of a substantial annulus of super-planetesimals in 
the absence of any remaining interior protoplanet which could 
accrete them. The stochastic fate of large individual bodies at 
late times can therefore overwrite and partially obscure system- 
atic trends in the data. In Scenario III, for example, if a single 
close encounter near the end of the simulation had resulted in a 
giant impact rather than a scattering into the star, a single inte- 
rior planet of (~ 12.5 M©) would have formed instead - another 
hot Neptune analogue. 

There are, however, many examples of hot Jupiters that 
have migrated further inward than 0. 1 AU. Orbits at ~ 0.05 AU 
are common and some objects have been found as close as 
~ 0.02 AU. It is less likely that interior planets would survive in 
such systems. If we assume that all remaining interior mass is 
accreted by the star (which is not always true, as demonstrated 
below) then the total surviving mass reduces to the exterior sur- 
viving mass which, as illustrated in Fig. 7, shows a clear trend 
with disk maturity. The later the migration episode the more 
disk mass will remain: from ~ 30% for our earliest scenario to 
~ 70% for our latest. None of our simulations support the con- 
jecture of a near-complete loss of solids from the swept zone. 



with disk maturity whereas that expelled into exterior orbits 
rises with disk maturity (see Fig. 7). The protoplanet mass frac- 
tion /pioto = 1 for the inner material is indicative of its rapid 
evolution, whereas the values of /proto ~ 0.6 - 0.9 for the outer 
material are as much influenced by the preferential tendency 
for protoplanets to be widely scattered as it is by their previous 
growth. The implication is that scattering is noticeably more 
effective in older, less dissipative disks. Giant planet migra- 
tion through a protoplanet/planetesimal disk rapidly advances 
its evolution (as measured by /proto), both by speeding up accre- 
tion and by fractionating objects according to the magnitude of 
their damping, but disks in an oligarchic stage of growth have 
more of a tendency to push material ahead of the giant and disks 
undergoing giant impact type growth allow a greater amount of 
mass to escape into external orbits. 

The amount of mass accreted by the giant or ejected are mi- 
nor in all cases but the picture is complicated by the statistics of 
mass accreted by the central star. In the last three scenarios, this 
is significant. However, the high values in Scenarios III & IV 
result mainly from the impact of a single massive protoplanet 
during the final energetic phase of accretion within 0.1 AU. 
The loss of mass to the star in Scenario V resulted from the 



3.3. The interior planets. 

In all scenarios, baiTing the latest, giant planet migration was 
found to stimulate accretion within the portion of the disk shep- 
herded inward. By the end of the simulations, this mass had 
accumulated into one or more massive planets within 0. 1 AU. 
Details of these planets and the giant, including their simula- 
tion ID, mass, a, e, and the presence of resonances are given in 
Table 4. 

A single interior planet was the most common result, but 
in one case there were three survivors. Their masses ranged be- 
tween ~ 2- 16 M© with semi-major axes between ~ 0.02 - 0.07 
AU. None of them remained in resonance with the giant even 
though some of their precursor bodies would have been origi- 
nally been pushed inward at first order resonances (see Fig. 4). 
These resonances were broken during the final accretion phase 
within 0.1 AU. In Scenario II however, the three planets that 
remain are all in resonant relation with each other The inner 
planet is in the 5:4 resonance with the middle planet and the 
5:3 resonance with the outer planet. The middle planet is in the 
4:3 resonance with the outer planet, giving a 5:4:3 commen- 
surability overall. It is possible that such a relationship could 
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Table 4. Interior surviving planets at 170000 years compared with the four currently known 'hot Neptune' systems. Data include the closest 
giant planet. 



Scenario 


ID 


Mass (Me) 


a(AU) 


e 


Resonances 


T 
1 




i J .UJ 


U.UJ J 




M 

IN one 




GIAl 


160.3 


0.100 


0.003 




TT 
11 


r»T TQ 


Z,.JJ 




n IAS 

U. l^o 


'>-A wjith TPPQ S-^ with 

J-H- WlLIl LK^iiy, J.J Willi 1^^12,1 




ICE9 


6.69 


0.054 


0.071 


with TCFl 




ICEl 


J.JO 




0.022 






GIAl 


160.9 


0.100 


0.005 




III 


ICEl 


4.89 


0.023 


0.110 


None 




GIAl 


160.3 


0.100 


0.001 




IV 


OLI28 


2.86 


0.036 


0.035 


None 




GIAl 


159.8 


0.100 


0.002 




GJ436 


b 


21.3 sin I 


0.028 


0.12 


None 


55 Cancri 


e 


14.3 sin! 


0.038 


0.174 


None 




b 


267.3 sin ; 


0.11 


0.02 




yti Arae 


d 


13.4 sin; 


0.09 


0.0 


None 




b 


541.0 sin;' 


1.5 


0.31 




GJ 876 


d 


5.89 sin/ 


0.021 


0.0 


None 




c 


177.9 sin/ 


0.13 


0.27 





-Scenarto I- 



how tt O Scervario II - 

-KH Q Scenario III - 
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Fig. 8. Comparison of computer generated interior planets with the 
interior regions of four known "hot Neptune" systems. 



act to stabilize the orbits of these planets, but they are closely 
spaced, the orbits of OLD and 1CE9 cross, and the giant acts 
to perturb the system, so we suspected accretion here to be in- 
complete. Running the inner system of Scenario II for an addi- 
tional 1 .0 Myr resulted in a prompt giant impact between the 
inner pair (1CE9 & OLI9, see Table 4), followed by a longer 
phase of interaction of the two 9.02 and 5.58 Me survivors. 
Their orbits gradually became more elliptical, especially that 
of the lighter outer planet (ICEl). Seven close encounters fol- 
lowed causing an outward scattering and a further excitation 
of the outer planet's eccentricity. Just before the million years 
was up, ICEl encountered the giant for the first time and 18 
close encounters later they collided. Thus, the final outcome 



for Scenario II was a single remaining inner planet, separated 
from the giant by 12.9 mutual Hill radii with a mass of 9.02 M^, 
fl = 0.0465 AU,e = 0.132. 

Generation of these massive interior planets by these sim- 
ulations is particularly interesting as three examples of short 
period Neptune-mass objects have been recently discovered in 
the systems GJ 436 (Butler et al. 2004), 55 Cancri (McArthur 
et al. 2004) and yu Arae (Santos et al. 2004) and a short pe- 
riod planet about half as massive may also have been detected 
in the GJ 876 system (Rivera et al. 2005). These systems are 
compared with those generated here in Table 4 and Fig. 8. In 
three of these natural systems the inner Neptune-mass planet is 
accompanied by more than one outer giant and, in the cases of 
55 Cancri and GJ 876, the innermost giant is placed close to 
where our simulated giant ends its migration. The results of the 
simulations do have a particular resemblance to reality in these 
two cases. The best matches are given by Scenario 1, where 
the mass and orbital radius of the interior planet are similar to 
the m sin / and a of 55 Cnc e, and Scenario III where the re- 
semblance is closer to the configuration of GJ 876 d & c. One 
might speculate therefore that, rather than hot Neptune type 
planets forming far out in the disk and migrating inward to their 
present location, they might, as illustrated here, form at these 
locations from disk material shepherded and compacted by a 
migrating giant. 

Against this proposition is the case of GJ 436 where the 
hot Neptune appears unaccompanied by an exterior giant; al- 
though Butler et al. (2004) did detect a linear velocity trend 
in their data implying the possible existence of a more distant 
companion. The primary of this system is a ~ 0.4 Mq red dwarf 
star which might affect the comparison. Giant planets appear 
to be rare in red dwarf systems and it may be that they do not 
form efficiently from the lower mass protoplanetary disks ex- 
pected around low mass stars (Laughlin et al. 2004). GJ 436 b 
may therefore be this system's largest planet, rather than a sec- 
ondary object, which could indeed have formed at large radius 
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Fig. 9. Scenario I: collision of the inner 
planet with the central star as the giant 
moves inward to 0.05 AU. Resonant angles 

for the 2:1 resonance are read on the left 
hand axis: light grey symbols plot the an- 
gle (p\ = 2/1' -A-ru' and dark grey symbols 
plot the angle 02 = 2A' -A-vj where A', 
and A, ru are the mean longitudes and lon- 
gitudes of periastron for the outer and inner 



17Sm irosu ITSW tntUB planet respectively. The semi major axis, pe- 

Vipni riastron and apastron for both the giant and 

the inner planet are plotted as black lines 
and are read on the right hand axis. 



before migrating in. However, not all red dwarf stars lack gi- 
ant planets. In the one known case where giants are present 
(GJ 876) and past migration may have played a role in their 
current configuration (e.g. Snellgrove et al. 2001), an interior 
half-hot Neptune appears to be present. The hot Neptune in the 
H Arae system is accompanied by a giant, but it is situated much 
further out (at 1 .5 AU) than the final location of the giant in the 
simulations. No simulations have yet been performed to evalu- 
ate the outcome of stopping the giant migration at 1.5 AU, but 
at the point in the runs where the giant passes through 1 .5 AU 
(at ~ 140 000 years) about 60% of the original mass of the disk 
is compacted within ~ 1 AU. This is enough mass to assem- 
ble a ^ 13 M© hot Neptune over a longer period although, 
since ju Arae d and n Arae b are separated by ~ 22 mutual 
Hill radii, one might expect some additional smaller planets to 
have formed and perhaps to have survived to the present. Less 
massive interior planets of ~ 2 - 7 are produced in the sim- 
ulations and may exist in nature as smaller versions of the hot 
Neptunes already discovered. GJ 876 d may represent the first 
discovery of a planet in this mass range. However, such objects 
would only produce a ~ 1 - 3 m s"' stellar radial velocity at 
sin ; = 1, a value that is close to the observation Umit, so others 
will have escaped detection so far. 

Another reason for the apparent rarity of hot Nepture type 
objects could be that many hot Jupiter systems are more com- 
pact than 0.1 AU; a typical example being 51 Pegasi b: m sinj 
0.5 Mj, a ^ 0.05 AU (Mayor & Queloz 1995). For a planet 
shepherded well within this distance, significant eccentricity 
excitation by the giant companion may cause it to impact the 
star. 

To examine this possibility, Scenarios I-IV were run for 
an extra 300 years, allowing the giant in each case to migrate 
further in to stop at 0.05 AU. In Scenarios I & II all the inte- 
rior planets were driven into the central star. The mechanism at 



work is illustrated for the case of Scenario I in Fig. 9 where the 
semi major axis, periastron and apastron distances, for both the 
giant and terrestrial planet, and the resonant angles for the 2:1 
resonance, are all plotted against time. The orbit of the interior 
planet is initially undisturbed, but when the giant has moved 
inward to 0.087 AU the planet is captured into the 2: 1 reso- 
nance as can be seen from the libration of the resonant angles. 
From this point the planet is pushed in ahead of the giant at 
the 2:1 resonance, its eccentricity increasing progressively. By 
the time the planet has reached a x 0.035 AU, its eccentricity 
has increased to e « 0.6 and impact with the star occurs at pe- 
riastron. The events in Scenario II were similar: the migrating 
giant in this case caused the three interior planets present to ac- 
crete each other, capturing the single 14.6 survivor into the 
2:1 resonance. From here, evolution proceeded as in Fig. 9, the 
planet eventually hitting the star through having been forced 
into a tighter, more elongated, orbit. 

In Scenarios III & IV, both the interior planets survived this 
additional migration by the giant. This is because their initial 
orbits were closer to the star (see Table 4) so the planet only 
becomes captured at the 2:1 resonance much later, or not at 
all. The interior planet in Scenario III remained stable as the 
2: 1 resonance did not reach its location and sweeping higher 
order resonances had no apparent effect. Survival of such a 
planet for the long term is therefore probable, since orbital de- 
cay due to tidal interaction with the star is expected to be very 
slow (see below). In Scenario IV, the interior planet was cap- 
tured into the 2: 1 resonance late: at the point where the giant 
reached 0.057 AU. Its orbit was compressed and elongated to 
a = 0.032 AU and e = 0.35 but the progressive increase in ec- 
centricity ceased when the giant finished migrating and impact 
with the star was avoided (see Fig. 10). 

The effect of tidal interactions with a slowly rotating star 
will cause the lower mass planet obtained in Scenarios III and 
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Fig. lO.Scenario IV: survival of the inner 
planet as the giant moves inward to 0.05 AU. 
Resonant angles for the 2:1 resonance are 
read on the left hand axis: hght grey sym- 
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IV to migrate inward slowly. For a circular orbit tidal dissipa- 
tion occurs within the star only, with turbulent convection in 
the stellar envelope being responsible for dissipating the tidally 
induced motions. The estimated orbital evolution time in this 
case is given by (Terquem et al. 1998): 



4 I Ma 
torb ^ 2.8 X 10-4 



1 day 



Gyr. 



(10) 



where P is the orbital period and lUp is the mass of the planet. 
For a 15 M® planet with an orbital period of 1.5 days this gives 
a orbital decay time of ~ 36 Gyr - i.e. comfortably longer than 
the age of the universe. For a planet on an eccentric orbit, tidal 
dissipation within the planet becomes important due to more 
effective dissipation within the solid body (e.g. Goldreich & 
Soter 1966). The raising and dissipation of tides within the 
planet leads to eccentricity damping, and the low moment of 
inertia of the planet ensures that it maintains near-synchronous 
rotation. This also means that effective removal of orbital angu- 
lar momentum can only be achieved through tidal dissipation 
in the star. We therefore expect that a planetary system consist- 
ing of an exterior gas giant planet and an interior rocky planet, 
orbiting in close proximity to a solar type star, will evolve such 
that the rocky planet orbit decays on a time scale on the order 
of that given by Eq. (10). The inner planet will have an eccen- 
tricity determined by a balance between tidal dissipation origi- 
nating in the planet itself and eccentricity excitation caused by 
the exterior giant. In general, the long term effects of this will 
be to make the orbits of both the inner terrestrial and outer giant 
planet more circular. This will be enhanced by the tidal interac- 
tion between the gas giant and the central star which will also 
tend to circularize its orbit (e.g. Rasio et al. 1996), though we 
note that this effect is not particularly relevant for our simula- 
tions as the giant planet is assumed to maintain a near-circular 
orbit through tidal interaction with the gas disk. For the spe- 



cific case of Scenario IV, where the inner planet is in 2: 1 reso- 
nance with the gas giant, tidal dissipation will cause the rocky 
planet orbit to circularize at a smaller semi major axis, probably 
removing the inner planet from the 2: 1 resonance in the pro- 
cess. Using Eq. (10) we estimate that the rocky planet formed 
in Scenario III will spiral into its host star on a time scale of 
~ 54 Gyr, and that formed in Scenario IV on a time scale of 
^ 700 Gyr, again both comfortably longer than the age of the 
universe. 

If the disk shepherding and compaction scenario advanced 
here has some validity, then it suggests the formation of hot 
Neptunes or lesser massive terrestrial planets as a by-product of 
giant planet migration. However, if the giant comes to rest too 
close to the star (< 0.05 AU) then the effects of mean motion 
resonances may cause interior planets to hit the star. The best 
prospect therefore of detecting close orbiting terrestrial type 
planets of > might be in systems with a circumscribing 
giant at a ~ 0.05 - 1.5 AU. 

3.4. The exterior scattered disk. 

As the giant migrated through the inner disk it scattered ~ 30- 
70% of the disk mass into external orbits (see Figs. 5 & 7 
and Table 3). In each case, a diffuse and dynamically excited 
external disk was generated, composed predominantly of pro- 
toplanet material (/proto ~ 0.7 - 0.9): individual protoplanets 
having a wide range of mass and orbital parameters. Data for 
the external protoplanets are presented in Table 5, giving their 
number, mean and maximum masses and orbital inclinations, 
and their mean, minimum and maximum semi major axes and 
eccentricities. It can be seen from the data that the number of 
external protoplanets reduced and their masses increased with 
disk maturity. Similarly, there is a tendency for protoplanets 
from a more mature disk to be more widely scattered. The for- 
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Table 5. External surviving protoplanets at 170000 years. 



Scenario 


N 


"Iproto (Me) 


Wmax (Me) 


a(AU) 


amin (AU) 


«max (AU) 


e 




^max 


i ° 


'max 


I 


32 


0.13 


0.86 


4.59 


0.63 


9.31 


0.48 


0.19 


0.81 


5.65 


25.18 


11 


32 


0.13 


0.86 


7.04 


1.31 


28.24 


0.53 


0.24 


0.86 


3.64 


9.35 


111 


25 


0.19 


0.67 


5.81 


1.84 


13.07 


0.47 


0.16 


0.87 


4.33 


14.90 


IV 


23 


0.42 


0.99 


8.18 


1.21 


33.5 


0.53 


0.13 


0.89 


2.71 


11.54 


V 


16 


0.55 


1.34 


10.26 


2.42 


40.97 


0.52 


0.10 


0.96 


4.62 


11.33 



mer trend is primarily due to previous accretion, whereas the 
latter is a result of the decreased dynamical friction operating 
at later times. 

The ejecta that comprised the scattered disk were spread 
over a much wider volume than their original location at 
< 4 AU. Planetesimal orbits were damped quite rapidly by gas 
drag to form a thin (~ 1 - 2 M^) external disk with a sur- 
face density ^ a few percent of the pre-existing Ss. Protoplanets 
were often in highly eccentric orbits, passing well beyond the 
confines of the original disk but with their periastra still located 
close to the location of their scattering within 4 AU. Mean or- 
bital inclinations were comparable to the solar system planets, 
but with a larger number of outhers as high as i » 25° (see 
Fig. 5). Thus, by selectively pushing planetesimals inward and 
widely scattering its external ejecta, the migrating giant par- 
tially evacuates a cavity within its swept zone. 

Further accretion in this disk will therefore be characterized 
by low Ss and high random velocities, reducing both the mass 
available and the effect of gravitational focussing. In some col- 
Usions, impact velocities could be high enough to cause disrup- 
tion of the protoplanets rather than accretion and a reversal of 
growth (Agnor & Asphaug 2004). Long evolution times are im- 
pUed for the mass contained in the scattered disk to rearrange 
itself into a smaller number of planets in stable orbits. This fi- 
nal configuration cannot be predicted from the juvenile stage 
illustrated in Fig.5. 

No original matter more distant than the giant's starting po- 
sition was modeled here. Objects in the external scattered disk 
traverse more widely than this and so could interact with matter 
in the outer disk beyond ~ 6 AU. A fresh supply of planetesi- 
mals would be encountered which could act to circularize and 
contract protoplanet orbits via dynamical friction (Thommes 
et al. 2002). Alternatively, encounters with other giant plan- 
ets remaining in the outer system could have a role to play in 
clearing material via ejection. The long term end product of the 
mass scattered by the migrating giant therefore depends par- 
tially on the nature of the outer disk: what other planets have 
formed there and its remaining population of small bodies. In 
a hot Jupiter system where no other gas giants have formed, 
we might speculate over the resulting planetary configuration 
at ~ 1 Gyr: a hot Neptune, or lesser massive terrestrial planet 
at ~ 0.05 AU, the giant at ~ 0. 1 AU, then from ~ 0.5 AU a suc- 
cession of Mars to Earth-mass planets, some still in eccentric 
or inclined orbits, extending as far out as ~ 5 AU. Beyond this, 
the content of the original outer disk would determine what is 
to be found. 

The ejection of low mass planets to large (~30 AU) semi- 
major axes may have some influence on the observed mor- 



Table 6. Total mass and number of protoplanets orbiting within, or 
crossing, the habitable zone (0.84 - 1.67 AU) at the end of the simu- 
lation. 



Scenario 


In Crossing Total 


I 


0.18 Me 2.13 M^ 2.31 Me 
N=l N = 8 N = 9 


n 


0.18 Me 1.79 Me 1.97 Me 
N=l N = 5 N = 6 


m 


1.28 Me 1.28 Me 

N=4 N=4 


IV 


1.17 Me 2.00 Me 3.17 Me 

N = 2 N = 3 N = 5 


V 


1.34 Me 1.34 Me 
N= 1 N= 1 



phology of the system in its debris disk phase. Systems such 
as Fomalhaut and e Eridani (Wyatt & Dent 2002; QuiUen & 
Thorndike 2002) have observational features that have been ex- 
plained by resonant trapping of planetesimals (Wyatt & Dent 

2002) or dust grains in mean motion resonances with a planet. 
We suspect that the dust trapping proposed by Quillen & 
Thorndike (2002) will also occur for lower mass planets (they 
considered nipianet ~ 30Me), but is likely to occur for reso- 
nances that lie closer to the planet than the 3:2 resonance that 
was dominant in their study. 

What of the probability of one of these surviving planets re- 
siding in the system's habitable zone? Habitable zones are dy- 
namically stable in hot Jupiter systems (Menou & Tabachnik 

2003) so a planet forming there would have a stable orbit if 
sufficiently well-spaced from neighbors. The simulations pre- 
sented here cannot provide any numerical estimate that might 
address this question: the simulations have not been run for 
long enough and, in any case, ignore the strong potential influ- 
ence of an outer disk. However, protoplanets are found within, 
or with their orbits passing though, the habitable zone (taken to 
be 0.84 - 1.67 AU; Kasting et al. 1993) at the end of each run 
(see Table 6). Thus, some mass is available for forming a com- 
pleted planet in the right place, especially in view of the fact 
that damping of the protoplanet's orbits from interaction with 
outer disk material, coUisional damping during ensuing accre- 
tion, destructive collisions, and gas drag on remaining planetes- 
imals and debris could act to return some material into closer 
orbits. Assuming eventual re-circularization of orbits with con- 
servation of angular momentum, the data in Table 6 convert 
to those illustrated in Fig. 1 1 which shows, for each scenario, 
the number of protoplanets and their total mass predicted in 
the habitable zone. Protoplanets occupy the habitable zone in 
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Fig. 11. Number and total mass of protoplanets predicted in the hab- 
itable zone assuming orbital re-circularization with conservation of 
angular momentum. 



three out of the five cases and, in two scenarios, more than an 
Earth-mass of material is present. Thus, since there is dynam- 
ical room available within the habitable zone of a hot Jupiter 
system, it is perhaps more likely than not that long term evo- 
lution of the scattered external disk could result in a terrestrial 
planet being located there. 

4. Caveats and future model improvements 

Inevitably a number of assumptions have been used in 
constructing the models presented in this paper, with the 
consequence that potentially important physical processes 
have been omitted. Here we discuss the possible implications 
of these for our results, whilst noting that work is underway to 
include them in future models: 

(0 Gap formation & cavity clearing: A giant planet is 
expected to form an annular gap in the gas disk centred 
around its orbit. As a Jupiter mass planet migrates inward, 
hydrodynamic simulations indicate that the inner disk becomes 
depleted of gas, forming a low density cavity there (e.g. 
Nelson et al. 2000). This arises in part because the viscous 
time scale in the inner disk is shorter than the migration time. 
The depletion is probably enhanced artificially by the use of 
an outflow boundary condition at the disk inner edge, which 
for computational reasons is located at a radius much larger 
than the expected surface of the central star. At this time it is 
not known accurately to what degree the inner disk is depleted. 
We have assumed an undepleted disk, and the effects of this 
will be the subject of future investigation. We note that our use 
of a 0.5 My planet will lead to significantly less gas depletion 
than obtained in hydrodynamic simulations performed using 
more massive planets. 

The effect of depleting the gas disk interior to the giant 
planet will be to reduce the gas drag experienced by the 
planetesimals. Reducing the dissipation level is likely to lead 



to greater scattering of bodies into the outer disk by the giant 
planet, which may reduce the efficiency of planet formation in 
the inner disk. We note, however, that the loss of planetesimals 
into the central star due to gas drag will also diminish. Models 
similar to those presented here, but coupled to an evolving gas 
disk model, will be the subject of a future publication. 
(//) Gas disk removal: Implicit within our models is a 
mechansim for preventing the gas giant planet from migrating 
all the way into the central star. While the presence of a 
magnetospheric cavity has been cited as a possible stopping 
mechanism for hot Jupiters (Lin et al. 1996), we note that 
planets are found to exist with a range of semi major axes for 
which this mechanism cannot be invoked. A more plausible 
reason for migration halting is the removal of the gas disk 
during migration. The removal of the outer gas disk would 
remove a source of dissipation for the scattered protoplanets 
and planetesimals, increasing the accretion time scale for any 
planets forming exterior to the giant. 

(///) Giant planet eccentricity evolution. Our model assumes 
that type II migration is associated with an eccentricity 
damping which maintains a near-circular orbit. However 
the question of whether a giant planet has its eccentricity 
damped or excited by the protoplanetary disk remains a point 
of ongoing debate (Papaloizou et al. 2001; Goldreich & Sari 
2003; Ogilvie & Lubow 2003). In the absence of eccentricity 
damping, interaction of the giant with the solids disk may 
cause modest excitation, but not at a level to significantly effect 
our results. If the gas disk was to drive significant eccentricity 
increase then this would probably lead to a stronger interaction 
between the giant and the solids disk resulting in greater 
scattering and a more efficient clearing of material. 

(iv) Planetesimals exterior to the gas giant: There will be a 
population of planetesimals exterior to the giant planet whose 
size distribution is unknown. Planetesimals with radii < 1 km 
will migrate in behind the gas giant due to gas drag, providing 
a source of material that may be accreted by outwardly scat- 
tered protoplanets that reside interior to 5 AU, and which may 
contribute to the damping of their inclination and eccentricity. 

(v) Planetesimal size evolution: The model we have used as- 
sumes that a population of 10 km sized planetesimals coexists 
with a population of larger protoplanets at the beginning of gas 
giant migration. Collisions between planetesimals may lead 
to both their growth and fragmentation, with a range of sizes 
developing. At present these processes are not included in our 
model, and their effect on the outcome of our calculations 
is unclear We simply note that smaller bodies produced by 
fragmentation will experience larger gas drag forces, leading 
to a more rapid in-spiral toward the central star, and reduced 
eccentricity that assists rapid accretion by the protoplanets. 
The build-up of larger bodies will lead to reduced effectiveness 
of gas drag, and a greater probability of scattering. 

ivi) Type I migration: We have neglected the effects of type 
I migration which operates for non-gap forming sub-Jovian 
planets (Ward 1997). This may become important for bodies 
more massive than ~ 1 M^, causing inward migration and 
strong eccentricity damping. For our model disk parameters 
type I migration proceeds faster than the gas giant migration 
for bodies more massive than ~ 5 M^. We note that such 
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bodies usually fomi during a rapid burst of accretion as the 
planet approaches the star, so their formation is unlikely to 
be greatly affected by type I migration. Their subsequent 
evolution will depend on the gas density in the inner disk, 
which at present is unknown, as discussed in point (/) above. 
We further comment that a consistent model of giant planet 
formation via the core-instabiUty model is difficult to achieve if 
type I migration operates as efficiently as current calculations 
suggest. Our model of a gas-giant forming and then slowly 
migrating inward implicitly contains the assumption that type 
I migration does not operate efficiently. We note that recent 
simulations of low mass planets in turbulent disks suggest 
that some low mass planets may avoid rapid inward migration 
(e.g. Nelson & Papaloizou 2004; Nelson 2005), so the role 
played by type I migration during planet formation is unclear 
at present. 

Type I disk-planet interactions are also thought to include 
significant eccentricity damping on a timescale ~100 times 
shorter than the migration timescale (Papaloizou & Larwood 
2000), an effect that could be relevant to the picture we have 
presented. The influence of this in simulations would be to 
exert additional dissipation on the protoplanets and hence 
increase their probability of being shepherded by the giant. 
Such behavior would have its greatest impact in more evolved 
systems where dissipation is weaker, such as our Scenarios IV 
and V. Type I damping would also effect the evolution of the 
scattered external disc and would be an additional influence 
acting to return some planet-forming material into closer 
orbits. Its integrated effect however would depend on how 
much of the nebula's gas remains and how much longer it lasts 
for. Eccentricity damping from a thin residual gas disk could 
act to circularize the orbits of any final planetary configuration 
(Kominami & Ida 2002), whereas if ample gas persists and 
type I migration is also in force, a pileup of mass could 
also occur exterior to the giant as migrating protoplanets are 
entrained in its exterior mean-motion resonances (Thommes 
2005). It is possible therefore that a modest degree of type I 
eccentricity damping could promote, rather than hinder, the 
accretion of a planetary system external to the orbit of the 
giant planet. 



5. Conclusions. 

Whilst this study has addressed the effect of giant planet mi- 
gration through inner disks of varying maturity, we have not 
examined the effect of changing some of the basic parameters 
that could influence accretion, shepherding and scattering be- 
havior. These might include varying the giant's mass, its migra- 
tion time (tv) and the mass of the nebula (/neb)- Such a study 
will be the subject of future work, which will also address some 
of the issues raised in Section 4 concerning extensions to the 
basic physical model. It is likely that varying the mass of the 
nebula would result in a change in accretion time scales and a 
corresponding increase or decrease in the mass of the surviving 
planets, but not in their elimination entirely. It would also affect 
the level of dissipation in the disk which influences the dynam- 
ical partitioning of disk mass brought about by the migration. 



Increasing the mass of the giant would raise the overall effect 
of scattering, as would extending the migration time which in- 
creases opportunities for repeated and incremental scatterings 
(Mandell & Sigurdsson 2003). It seems reasonable to speculate 
that varying any of these parameters within reasonable limits 
would still produce a result that, whilst varying in detail, re- 
mains consistent with our established picture. 

There are three principal conclusions that arise from this 
work. 

1 . Migration of a giant planet through an inner disk partitions 
the mass of that disk into internal and external remnants. 
The fraction of disk mass in either remnant is dependent on 
the level of dissipation and is thus sensitive to the maturity 
of the disk material at the time of the migration episode. 
Late migration favors the escape of more material into ex- 
ternal orbits. The survival of an inner remnant is also sensi- 
tive to the final position of the giant at the end of migration, 
this becoming increasingly unlikely in hot Jupiter systems 
with a ^ 0.05 AU. The concept that giant planet migration 
would eliminate all the mass in its swept zone is not sup- 
ported by our results. 

2. Hot Neptunes and lesser massive terrestrial planets are a 
possible by-product of type II migration, being formed 
from an inner system disk compacted by a migrating giant. 
Future searches of hot Jupiter systems for radial velocity 
signals close to the current detection limit might uncover 
more examples of these planets. 

3. Our results are supportive of the eventual accumulation of 
a number of terrestrial planets orbiting exterior to the gi- 
ant, including within the system's habitable zone. Thus, the 
early evolution and the final architecture of Hot Jupiter sys- 
tems does necessarily eUminate their possibility of hosting 
Earth-like planets. 
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